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Abstract. We give a detailed microlocal study of X-ray transforms over geodesics-like families of curves 
with conjugate points of fold type. We show that the normal operator is the sum of a pseudodifferential 
operator and a Fourier integral operator. We compute the principal symbol of both operators and the 
canonical relation associated to the Fourier integral operator. In two dimensions, for the geodesic transform, 
we show that there is always a cancellation of singularities to some order, and we give an example where 
that order is infinite; therefore the normal operator is not microlocally invertible in that case. In the case of 
three dimensions or higher if the canonical relation is a local canonical graph we show microlocal invertibility 
of the normal operator. Several examples are also studied. 



The purpose of this paper is to study X-ray type of transforms over geodesics-like families of curves with 
caustics (conjugate points). We concentrate on the most common type of caustics — those of fold type. Let 
70 be a fixed geodesic segment on a Riemannian manifold, and let / be a function which support does not 
contain the endpoints of 70. The question that we are trying to answer is the following: what information 
about the wave front set WF(/) of / can be obtained from the assumption that (possibly weighted) integrals 



of / along all geodesies 7 close enough to 70 vanish (or depend smoothly on 7)? Since X has a Schwartz 
kernel with singularities of conormal type, Xf could only provide information for WF(/) near the conormal 
bundle A/"* 70 of 70. If there are no conjugate points along 70, then we know that WF(/) n A/"* 70 = 0. 
This has been shown, among the other results, in [5J[T7] in this context. It also follows from the microlocal 
approach to Radon transforms initiated by Guillemin 5] when the Bolker condition (in our case that means 
no conjugate points) is satisfied. Then the localized normal operator N x := X*\X, where \ is a standard 
cut-off near 70 is a pseudo-differential operator (^DO), elliptic at conormal directions to 70. If there are 
conjugate points along 70, then N x is no longer a ^DO. One of the goals of this work is first to study the 
microlocal structure of N x in presence of fold conjugate points, and then use it to see what singularities can 
be recovered. That would also allow us to tell whether the problem of inverting X is Fredholm or not, and 
would help us to determine the size of the kernel, and to analyze the stability and the possible instability of 
this problem. 

Geodesic X-ray transforms have a long history, generalizing the Radon type X-ray transform in the 
Euclidean space, see, e.g., [7]. When the weight is constant, and (M,g) is a simple manifold with boundary, 
uniqueness and non-sharp stability estimates have been proven in [T31 Q31 using the energy method. 
Simple manifolds are compact manifolds diffeomorphic to a ball with convex boundary and no conjugate 
points. The uniqueness result has been extended to not necessarily convex manifolds under the no-conjugate 
points assumption in j4] . The authors used microlocal methods to prove a sharp stability estimate in [16] for 
simple manifolds and uniqueness and stability estimates for more general weighted geodesic-like transforms 
without conjugate points in [BJ. The X-ray transform over magnetic geodesies with the simplicity assumption 
was studied in [3] . Many of those and other works study integrals of tensors as well and the results for tensors 
of order two or higher are less complete. 
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1. Introduction 
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The authors considered in |17) the X-ray transform of functions and tensors on manifolds with possible 
conjugate points. Using the overdeterminacy of the problem in dimensions n > 3, we showed that if there 
exists a family of geodesies without conjugate points with a conormal bundle covering T* M, then we still 
have generic uniqueness and stability. In dimension two however that family has to be the set of all geodesies, 
and even in higher dimensions, |17j does not answer the question what is the contribution of the conjugate 
points to Xf. 



2. Formulation of the problem 

Let (M, g) be an n-dimensional Riemannian manifold. Let exp p (v), where (p,v) € TM, be a regular 
exponential map, see section [3j where we recall the definition given by Warner in |20j . The main example is 
the exponential map of g or that of another metric on M or other geodesic-like curves, for example magnetic 
geodesies, see also [5]. Let re be a smooth function on TM \ 0. We define the weighted X-ray transform Xf 

by 

(2.1) Xf(p,9)= f K(e Wp (t9),e±p p (t0))f(en> p (t9))dt, (p,0) € SM, 

where we used the notation 

exp(tv) = — exp(to). 

The t integral above is carried over the maximal interval, including t — 0, where exp(t9) is defined. The 
assumptions that we make below guarantee that this interval remains bounded. 

Let (po^o) € TM be such that v = vo is a critical point for exp po (t>) (that we call a conjugate vector) 
of fold type, see the definition below. Let go = ex P Po ( u o)- Then our goal is to study Xf for p close to po 
and 9 close to 9q : — ^o/l^ol under the assumption that the support of / is such that vq is the only conjugate 
vector v at po so that exp po (v) S supp/. Note that Vq can be written in two different ways as t9o, \9q\ = 1, 
with ±i > 0, and we chose the first one. The contribution of the second one can be easily derived from our 
results by replacing 9q by —9q. 

Instead of studying X directly, we study the operator 



(2.2) 



Nf(p)= / K t(p,6)Xf(p,e)da p (e) 

JS„M 



S V M 



K»(p,0)/e (exp p (td),exp p (t6)) f(exp p (t9))dtda p (9), 



with some smooth re" localized in a neighborhood of (po, 9q). Here da p {9) is the induced Riemannian surface 
measure on S p (M). When exp is the geodesic exponential map, there is a natural way to give a structure 
of a manifold to all non-trapping geodesies with a natural choice of a measure, see section [5j The operator 
X can be viewed as map from functions or distributions on M to functions or distributions on the geodesies 
manifold. Then one can define the adjoint X* with respect to that measure. Then the operator X*X is of 



the form (2.2) with re" = re, see (5.1). The condition that supp re" should be contained in a small enough 
neighborhood of {po,9 Q ) can be easily satisfied by localizing p near p , and choosing supp re to be near 
ilpo,s a , ip ,e )- In the case of general regular exponential maps N is not necessarily X*X. 



A direct calculation, see [IB] and Theorem 5.1 shows that the Schwartz kernel of X*X in the geodesic 
case (see also [B] for general families of curves), is singular at the diagonal, as can be expected, and that 
singularity defines a ^DO of order —1 similarly to the integral geometry problem for geodesies without 
conjugate points. We refer to section [5] for more details. Next, singularities away from the diagonal exist 
at pairs (p, q) so that q = exp p (v) for some v, and d„ exp p is not an isomorphism (p and q are conjugate 
points). The main goal of this paper is to study the contribution of those conjugate points to the structure 
of X*X and the consequences of that. We actually study a localized version of this; for a global version on 
a larger open set, under the assumption that all conjugate points are of fold type, one can use a partition of 
unity. 
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Let U be a small enough neighborhoods of (po,&o) m SM. Let U be a small neighborho od o f p so that 
U C n(U), where n is the natural projection on the base. Fix k" € Cq°(W). Let TV/ be as in (2.2 1, related to 
k", where k is a smooth weight. We will apply X to functions / supported in an open set V 9 po satisfying 
the conjugacy assumption of the theorem below, see Figure [I] Our goal is to study the contribution of a 
single fold type of singularity. Let S C M x M be the conjugate locus in a neighborhood of (po,qo), see 
section [3] Finally, let 70 = 7 Po ,0 o (i), t & I, be the geodesic through {po,@o) defined in the interval / 3 0, 
with endpoints outside V. 




Figure 1. 

The first main result of this paper is the following. 

Theorem 2.1. Let Vq = \vq\9q be a fold conjugate vector at po, and let N be as in fi2.0jf . Let vq be the only 
singularity o/exp po (u) on the ray {exjp p (t9o), t G /} (1 V. Then if U (and therefore, U) is small enough, 
the operator 

N : C™(V) — > C^°(U) 

admits the decomposition 

(2.3) N = A + F, 

where A is a ^ DO of order —1 with principal symbol 



(2.4) 



a p (A)(x, = 2tt / 6(t(e)) («*«)(«, 0) da x (6), 



and F is an FIO of order —n/2 associated to the Lagrangian 7V*E. In particular, the canonical relation C 
of F in local coordinates is given by 

(2.5) C = {(p,£,q,ri), (p,q) G S, £ = —r]idexpp(v)/dp,r) G Cokerd^ exp p (v), det d„ exp )9 (w) = 0} . 

7/exp is £/ie exponential map of g, then C can also be characterized as Af*Y,' , where TVS is as in j[4-17 ), 
and the prime means that we replace r\ by —r\. 

It is easy to check that C above is invariantly defined. 

In section [9] we show that in dimension 3 or higher in the case that C is a local canonical graph the 
operator N is microlocal invertible. In two dimensions, in the geodesic case, we show that there is always a 
loss of some derivatives at least when the curves are geodesies. We study in detail the case of the circular 



Radon transform in two dimensions in section 10 and show that then N is not microlocally invertible. 



3. Regular exponential maps and their generic singularities 

3.1. Regular exponential maps. Let M be a fixed n-dimensional manifold. We will recall the definition 
of Warner [20 of a regular exponential map at p G M. We think of it as a generalization of the exponential 
map on a Riemannian manifold, by requiring only those properties that are really necessary for what follows. 
For that reason, we use the notation exp p (y). In addition to [2U] , we will require exp p (-u) to be smooth in p 
as well. Let N p (v) C T v T p M denote the kernel of dexp p . Unless specifically indicated, d is the differential 
w.r.t. v. The radial tangent space at v will be denoted by r v . It can be identified with {sv, s G R}, where 
v is considered as an element of T v T p M. 
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Definition 3.1. A map exp p (v) that for eachp £ M maps v 3 T p M into M is called a regular exponential 
map, if 

(Rl) exp is smooth in both variables, except possibly at v = 0. Next, dexp p (it>)/di ^ 0, when v =/= 0. 
(R2) The Hessian d 2 exp p (v) maps isomorphically r v x N p (v) onto T oxp ( V ^M /dexp p (T v T p M) for any 

n/Oin T p M for which exp p (t>) is defined. 
(R3) For each v € T p M \ 0, there is a convex neighborhood U of v such that the number of singularities 

of exp p , counted with multiplicities, on the ray tv, t € R in U, for each such ray that intersects U, 

is constant and equal to the order of v as a singularity of exp p . 

An example is the exponential map on a Riemannian (or more generally on a Finsler manifold), see |20) . 
Then (Rl) is clearly true. Next, (R2) follows from the following well known property. Fix p and a geodesic 
through it. Consider all Jacobi fields vanishing at p. Then at any q on that geodesic, the values of those 
Jacobi fields that do not vanish at q and the covariant derivatives of those that vanish at q span T q M . 
Also, those two spaces are orthogonal. Finally, (R3) represents the well known continuity property of the 
conjugate points, counted with their multiplicities that follows from the Morse Index Theorem (see, e.g., 
Thm 4.3.2]). 

We would need also an assumption about the behavior of the exponential map at v = 0. 
(R4) exp p (it>) is smooth in p, t, v for all p € M, \t\ <C 1, and n^O. Moreover, 

exp p (0) = p, and — exp p (tv) = v for t = 0. 

Given a regular exponential map, we define the "geodesic" 7 P .„(i), v ^ 0, by "y P ,v(t) = ex.p p (tv). We will 
often use the notation 

(3.1) q = cxp p (v) = 7 P ,„(1), w = -exp p (v) := -jp, v (l), 6 = v/\v\. 

Note that the "geodesic flow" does not necessarily obey the group property. We will assume that 



(R5) For q, w as in (3.1 ), we have exp„(«;) = p, exp q (w) = —v. 
This shows that in particular, (p, v) M- (q, w) is a diffeomorphism. If exp is the exponential map of a 
Riemannian metric, then (R5) is automatically true and that map is actually a symplectomorphism (on 
T*M). 

Remark 3.1. In case of magnetic geodesies, or more general Hamiltonian flows, (R5) is equivalent to time 
reversibility of the "geodesies." This is not true in general. On the other hand, one can define the reverse 
exponential map exp~(u>) = 7g,-u>(~l) in that case, see e.g. [3], near (qo,wo), and replace exp by exp~ in 
that neighborhood. Then (R5) would hold. In other words, (R5) really says that (p, v) i-> (q, w) is assumed 
to be a local diffeomorphism with an inverse satisfying (Rl) - (R4). 

3.2. Generic properties of the conjugate locus. We recall here the main result by Warner [20] about 
the regular points of the conjugate locus of a fixed point p. The tangent conjugate locus S(p) of p is 
the set of all vectors v € T p M so that dexp p (w) (the differential of exp p (u) w.r.t. v) is not an isomorphism. 
We call such vectors conjugate vectors at p (called conjugate points in [10]). The kernel of dexp p (v) is 
denoted by N p (v). It is part of T v T p M that we identify with T p M. In the Riemannian case, by the Gauss 
lemma, N p (v) is orthogonal to v. In the general case, by (Rl), it is always transversal to v. The images of 
the conjugate vectors under the exponential map exp p will be called conjugate points to p. The image of 
S(p) under the exponential map exp p will be denoted by £(p) and called the conjugate locus of p. Note 
that S(p) C T p M, while S(p) C M. We always work with p near a fixed Pq and with v near a fixed i>o- Set 
go = ex P Po ( w o)- Then we are interested in S(p) restricted to a small neighborhood of vq, and in S(p) near 
<7o- Note that S(p) may not contain all points near go conjugate to p along some "geodesic"; and may not 
contain even all of those along exp po (too) if the later self-intersects — it contains only those that are of the 
form exp p (w) with v close enough to v . 

Normally, dexp p (v) stands for the differential of exp p (w) w.r.t. v. When we need to take the differential 
w.r.t. p, we will use the notation d p for it, We write d v for the differential w.r.t. v, when we want to distinguish 
between the two. 
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We denote by E the set of all conjugate pairs (p, q) localized as above. In other words, E = {(p, q); q E 
E(p)}, where p runs over a small neighborhood of pq. Also, we denote by S the set (p,v), where v £ S(p). 

A regular conjugate vector v is defined by the requirement that there exists a neighborhood of v, 
so that any radial ray of T p M contains at most one conjugate point there. The regular conjugate locus 
then is an everywhere dense open subset of the conjugate locus that has a natural structure of an (n — 1)- 
dimensional manifold. The order of a conjugate vector as a singularity of exp p (the dimension of the kernel 
of the differential) is called an order of the conjugate vector. 

In [501 Thm 3.1], Warner characterized the conjugate vectors at a fixed p of order at least 2, and some of 
those of order 1, as described below. Note that in B\, one needs to postulate that N po (v) remains tangent 
to S(po) at points v close to vq as the latter is not guaranteed by just assuming that it holds at i>o only. 

(F) Fold conjugate vectors: Let vo be a regular conjugate vector at po, and let N Pq (vq) be one- 
dimensional and transversal to S(po). Such singularities are known as fold singularities. Then one 
can find local coordinates £ near vq and y near q so that in those coordinates, exp po is given by 

(3-2) y' = e, y n = (Cf- 

Then 

(3.3) S( Po ) = {t = 0}, N po (vo) = sp*a{d/d€ n }, E(po) = {y n = 0}. 

Since the fold condition is stable under small C°° perturbations, as follows directly from the defini- 
tion, those properties are preserved under a small perturbation of pg. 
(Bi) Blowdown of order 1: Let Do be a regular conjugate vector at po and let N Po (v) be one-dimen- 
sional. Assume also that N po (v) is tangent to S(po) for all regular conjugate v near vo. We call such 
singularities blowdown of order 1. Then locally, exp po is represented in suitable coordinates by 

(3.4) y' = Z', y n =eC- 
Then 

(3.5) S(po) = U 1 = 0}, A^>o)= span {d/dri, £(po) = {V = 2/" = 0}. 

Even though we postulated that the tangency condition is stable under perturbations of Vq, it is not 
stable under a small perturbation of po, and the type of the singularity may change then. In some 
symmetric cases, one can check directly that the type is locally preserved. 

(B/j) Blowdown of higher order: Those are regular conjugate vectors in the case where N Po (vq) is 
fc-dimensional, with 2 < k < n — 1. Then in some coordinates, exp po is represented as 



(3-6) 



(3-7) 



Then 



y % = C, i = l,...,n-k 
y 1 = ^C, i = n - k+ 1, . . . ,7i. 

S(Po) - U 1 = 0}, N pa (v ) = span {d/de~ k+ \ . . . , d/dC} 

e(po) = iv 1 = y n ~ k+1 = ■ ■ ■ = y n = o}. 



In particular, N Po (vq) is tangent to S(po). This singularity is unstable under perturbations of po, as 
well. A typical example are the antipodal points on S n , n > 3; then k = n — 1. 

The purpose of this paper is to study the effect of fold conjugate points to X. 



4. Geometry of the fold conjugate locus 

In this section, we study the geometry of the tangent conjugate locus S(p), and S respectively; and the 
conjugate locus E(p) and E, respectively. Recall that we work locally, and everywhere below, even if not 
stated explicitly, {p, v) belongs to a small enough neighborhood of (po, vq); (q, v) is near (go, wo)- We assume 
throughout the section that vq is conjugate vector at po of fold type. We also fix a non-zero covector rjQ at qo 
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as in (2.5 1, and let £ be the corresponding £ as in (2.5). We will see later that ^ 0. We refer to Figure [2] 
where w is not shown, and the zero subscripts are omitted. 
We start with properties of S(j>) and S. 

Lemma 4.1. 

(a) Let v £ S(p) be a fold conjugate vector. Then near q = cxp p (v) 7 is a smooth surface of 
codimension one, tangent to w :— — 7 p .„(l). 

(b) S is a smooth (2n — I) -dimensional surface in TM that can be considered as the bundle {S{p),p £ M} 
with fibers S (p) . 




Figure 2. A typical fold conjugate locus 



Proof. Consider (a) first. The representation (3.2 1 implies that locally, = exp p (5(p)) is a smooth 

surface of codimension one (given by y n = 0). Next, for v € S(p), the differential dexp p sends any vector 
to a vector tangent to S(p), as it follows from (3.2 1 again. In particular, this is true for the radial vector v 
(considered as a vector in T v T p M). This proves that w is tangent to E(p). 

The statement (b) follows from the fact that S is defined by detdexp p (u) = 0, and that detdexp p (i>) has 
a non-vanishing differential w.r.t. v. □ 

Remark 4.1. It is easy to show that in (a), 7 P)tI is tangent to of order 1 only. 

We define "Jacobi fields" along 7 Pj „ vanishing at p as follows. For any a £ T v T p M, set 

d 

J(t) = d[exp p (tv)](a) = a k —^e-xp^tv). 

Then J(0) = 0, J(0) = a, where J(T) = dJ(t)/dt. If J(l) = 0, then a direct computation shows that 

(4.1) j(l) =d 2 ex Pj »(ax v). 

When exp is the exponential map of a Riemannian metric, it is natural to work with the covariant 
derivative D t J(t) =: J'(l) instead of J(t). While they are different in general, they coincide at points where 
J {ft) = 0. 

The next lemma shows that the fold/blowdown conditions are symmetric w.r.t. p and q. 

Lemma 4.2. The vector vq is a conjugate vector at p$ of fold type, if and only if wq is a conjugate vector 
a t Qa of fold type. 
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Proof. Set wo — — j p0tVo (l), as in (3.1). Then p = exp 9o (uio). Assume now that a £ N po (vo). In some local 
coordinates, differentiate p — exp a (w) w.r.t. v in the direction of a; here q, w are viewed as functions of p, 



v. Then, using the Jacobi field notation introduced above in (4.1), we get 



because 



= dexp 9o Oo) ^a fe J^(p ,t'o)^ = dexp 9o (w ) J(l) 



k 9w , . k d d • 
01 a~k(P°> v °) = a «~T17 exp (to)(po,w ) = J(l). 
ov K ov h at t=i F 



By (R2), J(l) 0, so in particular, this shows that wq is conjugate at go, and J(l) € N qo {wo). Moreover, 
by (R2), the linear map 

(4.2) N p (v) 3 a = j(0) i-> j(l) := £ € J(0) = ,7(1) = 



defines an isomorphism between N p (v) and N q (w). Then (4.2 ) shows that wo is conjugate at go of multiplicity 
one. By (R3), applied to wq, it is also regular. 

We will prove now that wq is of fold type. Since it is regular and of multiplicity one, S(qo) near wq is 



a smooth (n — 1) dimensional surface either of type F, as in (3.3) or of type B±, as in (3.5). Assume the 
latter case first, then E(qo) is of codimension two, as follows from (3.5). In particular, using the normal 
form ( |3.4[ ), we see that in this case, one can find a non- trivial one-parameter family of vectors w(s) so that 
w(0) = wo and exp (Jo (w(s)) = pq. Then the corresponding tangent vectors at po would form a non-trivial 
one-parameter family of vectors v(s) so that exp po (w(s)) = go- That cannot happen, if vq is of type F, see 
(3.2), since the equation exp po (w) = go has (near Vq) at most two solutions. □ 



For (p, v) £ S, let a = ct(p, v) £ N p (v) be a unit vector. To fix the direction, assume that the derivative 
of det dexp p (w) in the direction of a, for v a conjugate vector, is positive. Here we identify in T v T p M and 
T p M. In the fold case, N p (v) is clearly a smooth vector bundle on TM near (po,Vo), and a is a smooth 
vector field. 

Lemma 4.3. For any fixed p near p , the map 

(4.3) S(p) 3«4 a(p, v) £ N p {v) 
is a local diffeomorphism, smoothly depending on p if and only if 

(4.4) d 2 exp po (v ) (N po (v Q ) \ x • ) \ TvoS(po) is of full rank. 

Proof. In local coordinates, we want to find a condition so that the equation 

a z d v i exp p (ii) = 

can be solved for v so that v = Vq for (p, a) = (po,cto), where ao = a(po,vo). Then v would automatically 
be in S(p). By the implicit function theorem, this is equivalent to 

det (d v a l d v i exp po (u)) =^ at v = vq. 

Choose a coordinate system near w so that d/dv n spans N Po (vq), and {d/dv 1 , . . . d/dv^ 1 } span T Vo S(po). 
Denote F(v) = exp po (v) and denote by Fi, Fij the corresponding partial derivatives. Greek indices below 
run from 1 to n — 1 . We have 

(4.5) d n F(v ) = 0, because d/dv n £ N Po (v ), 

(4.6) d a det(dF)(v ) = 0, because d/dv a is tangent to S(p ) at v , 

(4.7) d n det(dF)(v ) 7^ 0, by the fold condition, 

(4.8) c a d a F(v ) ^ 0, Vc^O, because c a d/dv a £ N Po (v Q ). 



We want to prove that det(d n dF)(vo) ^ if and only if (4.4) holds. That determinant equals 



(4.9) det{F ln ,F 2ni ...,F nn )(v ). 
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Perform the differentiation in (|4.6|. By (4.5), (4.8), 



det(F 1 ,...,F n _x,F na )(v ) = 0, \fa =4> F„a(uo) G span(.Fi(uo), . . . ,F n -i(v )). 



Similarly, (4.7) implies 



(4.10) det(F 1> ...,* , n _ 1 ,F nn )(uo)^0 => / F nn (v ) $ span^^o), ... ,F n _ x (v )). 



Those two relations show that (4.9) vanishes if and only if (F n i(vo), ■ ■ • i*n )n _i(i>o)) form a linearly dependent 



system, that is equivalent to (4.4). □ 



We study the structure of the conjugate locus E(p), E(g) and E next. Recall again that we work locally 
near po, vq and qo. 

Theorem 4.1. Let vq be a fold conjugate vector atpQ. 

(a) Then for any p near pq, E(p) is a smooth hypersurface of dimension n — 1 smoothly depending on p. 
Moreover for any q — exp p (v) £ E(p), T q M is a direct sum of the linearly independent spaces 

(4.11) T q M = T q E(p)®N q (w), 

and 

T q Yi(p) = Im dexp p (w), iV*E(p) = Cokerd„ exp p (v). 

Next, those statements remain true with p and q swapped. 

(b) E is a smooth (2n — 1) -dimensional hypersurface in M x M near (po,qo), that is also a fiber bundle 
E = {E(p), p £ M} with fibers E(p) (and also E = {E(g), q G Moreover, the conormal bundle 7V*S 
zs given by 

N*Z={(p,q,€,v); (p, q) £^, £,= Vtdexp p (v)/dp, rj £ Cokerd v exp p (v) 
where v = exp~ 1 (q) with exp p restricted to S(p) }. 



Proof. We start with (a). By the normal form ( |3.2[ ), also clear from the fold condition, the image of S(p) 
under dexp p (u) coincides with T q Y,(p). In particular, dexp p (u), restricted to S(p) is a diffeomorphism to its 



image. Relation (4.11) follows from (4.2) and (R2) 



Consider (b). We have (p, q) £ E if and only if there exists v (near vq) so that 
(4.13) q = exp p (v), det d„ exp p (?;) = 0. 

In some local coordinates, we view this as n + 1 equations for the 3n-dimensional variable (p, q, v) near 
(POi <70j v o)- We show first that the solution that we denote by L, is a (2n — 1) -dimensional submanifold. To 
this end, we need to show that the following differential has rank n + 1 at (pq, qo, Vq): 



(4.14) 



d p exp p (u) -Id d^expp^) 
v d p det d v exp p (v) d„ det d^ exp p (u) / 

The elements of the first "row" are n x n matrices, while the second row consists of three n- vectors. That the 
rank of the differential above is full follows from the fact that d v detdt, exp p (z;) ^ at (pq, vq), guaranteed 
by the fold condition. 

Set ir(p,q,v) = (p,q). We show next that tt(L) is a (2n — l)-dimensional submanifold, too. To this end, 
we need to show that d7r is injective on TL. The tangent space to L is given by the orthogonal complement 



to the rows of (4.14). Let us denote any vector in TL by p = (p p , p q , p v ). Then d7r(p) = (p p ,p q ). Our goal 



is therefore to show that p p = p q = implies p v = 0. Then (0,0, p v ) is orthogonal to the rows of (4.14), 
therefore, 

p\d v i exp p (v) =0, k = 1, . . . , n, p l v d v , det d„ exp p (v) = 0. 

The latter identity shows that p v £ N p (v), while the first one shows that p v £ Kerd„ exp p (?;). By the fold 
condition, p v = 0. 

This analysis also shows that the covectors v orthogonal to E are of the form v — (v p , v q ) with the property 



that (y p , v qi 0) is conormal to L. Since the conormals to L are spanned by the rows of (4.14), in order to get 
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the third component to vanish, we have to take a linear combination with coefficients a,, i 
so that 



. n and b 



(4.15) 



dq l 



ddet d„ exp„(v) 
b " . pK ' = 0, 



where q = exp p (v). Let 0^a£ N p (v). Multiply by a? and sum over j above to get that the ^-derivative of 
6det d„ exp p (w) in the direction of N p (v) vanishes. According to the fold assumption, this is only possible if 
b = 0. Then we get that a £ Cokerdt, exp p (v). Therefore the normal covectors to £ are of the form 

dq 1 



(4.16) 



■dpi 



a £ Coker d„ exp p (v), 



that proves (4.12). 



□ 



Theorem 4.2. Let vq be a fold conjugate vector at po. Let exp p be the exponential map of a Riemannian 
metric. 



(a) Then the sum in (4-H) is an orthogonal one, i.e., 

N q X(p) = N q (w). 



(b) Next, (4- 11) also admits the representation 
(4.17) 



A/"S ={(p, q, ct, f3); (p,q) = J'(0),/3 = — J'(l), where J is any Jacobi field 

along the locally unique geodesic connecting p and q with J(0) = J(l) = 0}. 

(c) is a graph of a smooth map (p,a) i-> (q,(3) if and only if condition (4-4) * s fulfilled. Then that 
maps is a local diffeomorphism. 

Remark 4.2. Note that for (p,q) <E S, the geodesic connecting p and q is unique, as follows from the 
normal form (3.2), only among the geodesies with 7(0) close to vq. Also, J is determined uniquely up to a 
multiplicative constant. Next, once we prove that E is smooth, then a € N p (v) and f3 £ N q (w) by (a) (see 
also (3.2)), but (4.17) gives something more than that — it restricts (a,j3) to an one-dimensional space. 

Remark 4.3. It is a natural question whether |J'(0)| = | J'(l)|. One can show that generically, this is not 
the case. 

Proof. By |12l Lemma IX. 3. 5], the conjugate of dexp p (w), w.r.t. the metric form is given by 

(4.18) (dcxp p (v))* = dexp q (w), 

where we use the notation (3.1 1. The normal to at q is in the orthogonal complement to the image of 
dexp p (v), that by (4.18) is Ker d exp^ (w) = N q (w). This proves (a). 

Then we get by (4.18 ), (4.15) (where 6 = 0) that a £ N q (w), where we identify the covector a with a 
vector by the metric. 

We will use now p21 Lemma IX. 3. 4]: for any two Jacobi fields J±, J2 along a fixed geodesic, the Wronskian 
(J[, J 2) — (Ji, J' 2 ) is constant. Along the geodesic connecting p and q, in fixed coordinates near p, let J be 
determined by J(0) = Cj, J'(0) = 0. Here ej has components Sj. If p and q are conjugate to each other, then 
J(l) is the equal to the variation dq/dpP , and this is independent on the choice of the local coordinates, as 
long as ej is considered as a fixed vector at p. Define another Jacobi field by J(l) = 0, J'(l) = a, where a 
is as in (4.16) but considered as a vector. Denote the field in the brackets in (4.16) by Xj. Then 

(4.19) A J = ( a ,J(l)) 

= <J'(1),J(1)> 

= <J'(1),J(1))-(J(1),J'(1)) 
= (J'(0),J(0))-(J(0),J'(0)) 
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This proves (4.17). 



The proof of (c) follows directly from Lemma 4.3 □ 



5. The Schwartz kernel of N near the diagonal and mapping properties of X and N 

5.1. The geodesic case. Let exp be the exponential map of the metric g. Then X is the weighted geodesic 
ray transform. One way to parametrize the geodesies is the following. Let H be any orientable hypersurface 
with the property that it intersects transversally, at one point only, any geodesic in issued from a point in 
14. For our local analysis, H can be an arbitrarily small surface intersecting transversally 7 Po ,u , so let us fix 
that choice. Let dVol# be the induced measure in H, and let v be a smooth unit normal vector field on H 
consistent with the orientation of H. Let T-L consists of all (p, 0) € SM with the property that p £ H and is 
not tangent to H, and positively oriented, i.e., (v, 9) > 0. Introduce the measure dp, = (n, 9) d Vol_y(p) der p (0) 
on %. Then one can parametrize all geodesies intersecting H transversally by their intersection p with H 
and the corresponding direction, i.e., by elements in Ti. An important property of dp is that it introduces a 
measure on that geodesies set that is invariant under a different choice of H by the Liouville Theorem, see 
e.g., PS]. 

The weighted geodesic transform X can be defined as in (2.1 ) for (p, 9) g H instead of (p, 9) GU because 
transporting (p, v) along the geodesic flow does not change the integral. Since we assumed originally that k 
is localized near a small enough neighborhood of 7 Po ,„ , we get that n is supported in a small neighborhood 
of (po, 0o ) m H- We view X as the following map 

X : L 2 (M) -> L 2 ("H, d/i), 

restricted to a neighborhood of (po, 0o)- This map is bounded, see |15j . and this also follows from our analysis 
of N. By the proof of Proposition 1 in [TB], X*X is given by 

(5.1) X*Xf(p) = — L= f f R(p,9) K (exp p (t9),eip p (t9))f(exp p (t9))dtda p (9). 

V det g{p) J s p m J 

We therefore proved the following. 

Prop osition 5.1. Let ex p be the geodesic exponential map. Let X be the weighted geodesic ray transform 
(2.1), and let N be as in \2.%fy , depending on k$ . Then 

X*X = N with «;» = it. 



(5.2) X*Xf{p) = -f=== I W(p,v)f(exp p (v))dVo\(v), 



Split the t integral in ( 5.1 ) in two: for t > and for t < 0, and make a change of variables (t, 9) n- (— t, 
in the second one to get 

1 

y/detgip) Jt p m 

where 

W = \v\- n+1 (R(p, v/\v\)k( exp («), exp Jv)/\v\) 
(5.3) V 

+ «(p, ~v/\v\)k( exp p (v), -exp p (v)/\v\) J . 

Note that |exp (t;)| = \v\ in this case. 



Next we recall a result in [TB]. Part (a) is based on formula (5.2) after a change of variables. 

Theorem 5.1 (|16j). Let exp be the exponential map of M. Assume that exp p : exp~ 1 (Af) — > M is a 
diffeomorphism for p near po . 

(a) Then for p in the same neighborhood ofpg, 



(5-4) X*Xf{p) = _ _ / A(p,q) 



1 r A ,„^ m det 5 ^/ 2 ) 



dpdq 



dq, 



where 

A(p, q) = R(p, - grad p)n(q, grad p) + k(j>, grad p)n{q, - grad p). 
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(b) X* X is a classical ^DO of order — 1 with principal symbol 



(5.5) 



a p (X*X)(x,0 = 2tt / S(aO))\<^,O)\ 2 da x (0), 



S. r M 



where £(0) = ^6 J , and 6 is the Dirac delta function. 



Note that the integral (5.4) is not written in an invariant form but one can easily check that writing it 



w.r.t. the volume form, the kernel is invariant. We also note that in the proof of Theorem 2.1 we apply the 



theorem above by restricting supp / and the region where we study N f to a small enough neighborhood of 
Po, where we there will be no conjugate points. This gives the ^DO part A of N in Theorem |2.1| 



Mapping properties of X. Let (x',x n ) be semigeodesic coordinates on H near xq. Then param- 
eterize the vectors near (xo,9q). We define the Sobolev space of functions constant along the flow, 
supported near the flow-out of (xq,9o) as the H s norm in those coordinates w.r.t. the measure dfi. We 
can chose another such surface H near q with some fixed coordinates on it; the resulting norm will the be 
equivalent to that on %. 

Proposition 5.2. With the notation and the assumptions above, for any s > 0, the operators 

(5.6) X : H ( ](V) — ► H S+1 / 2 {H), 

(5.7) X*X : H°{V) — -> H S+1 (V) 
are bounded. 

Proof. Recall first that the weight k localizes in a small neighborhood of (70,70)- Let first / has small 
enough support in a set that we will call Mq. Then Mq will be a simple manifold if small enough. Then we 
can replace H by another surface Ho that lies in Mo, and denote by Ho the corresponding H . This changes 
the original parameterization to a new one, that will give us an equivalent norm. 
Then, if s is a half-integer, 



l^/lliJa+l/2(K ) - ^ 



H<2s + 1 



(d% 4 ,xf,xf) 



= c 



E 

|a|<2s+l 



The term 9", ^,Xf is a sum of weighted ray transforms of derivatives of / up to order \a\. Then X*d", ^,X 
is a \I/DO of order \a\ — 1 because Mo is a simple manifold. That easily implies 

ll^/llff-+v. (WG ) < C\\f\\ H s. 

The case of general s > follows by interpolation, see, e.g., [151 Sec 4.2]. 

To finish a proof, we cover 70 with open sets so that the closure of each one is a simple manifold. Choose 
a finite subset and a partition of unity 1 = ^ZXj related to that. Then we apply the estimate above to 
each X\jf on the corresponding T-Lj. We then have finitely many Sobolev norms that are equivalent, and 



in particular equivalent to the one on T-L. This proves (5.6) 



To prove the continuity of X*X, we need to estimate the derivatives of X* X. We have that d a X*Xf is 
sum of operators X Ka of the same kind but with possibly different weights applied to derivatives of X f up 

0. For /, h in Cg>(V), |/ 



to order a , see (5.1). Let first s 



L 2 (V) 



= 1, we have 



In the last inequality, we used (5.6) that we proved already. This proves i5.7l for s = 



For s > 1, integer, we can "commute" the derivative in d a X*X with X*X by writing it as a finite sum 
of operators of the type X'^XpPpf, \/3\ < \a\, where Pp are differential operators of order f3. To this end, 



we first "commute" it with X* , as above, and then with X. Then we apply (5.7) with s — 0. The case of 
general s > follows by interpolation. □ 



12 



P. STEFANOV AND G. UHLMANN 



Remark 5.1. We did not use the fold condition here. In fact, Proposition |5.2| holds without any assumptions 
on the type of the conjugate points, as long as V is contained in a small enough neighborhood of a fixed 
geodesic segment that extends to a larger one with both endpoints outside V. Note that proving the 
mapping properties of X*X based on its FIO characterization is not straightforward, and we would get the 
same conclusion under some assumptions only, for example that the canonical relation is a canonical graph; 
that is not always true. 

Remark 5.2. A global version of Proposition |5 . 2| can easily be derived by a partition of unity in the phase 
space. Let (M, g) be a compact non-trapping Riemannian manifold with boundary. Let Mi be another such 
manifold which interior includes M, and assume that dMi is strictly convex. Such Mi always exists if dM 
is strictly convex. Let d-SMi denote the vectors with base point on dM pointing into M\. Then we can 
parameterize all (directed) geodesies with points in d-SM\, that plays the role of H above. Then for s > 0, 

X : H$(M) — > ff s+1 / 2 (9_S'Mi), X* X : H$(M) — ► ff s+1 (Mi) 

are bounded. 

5.2. General regular exponential maps. Let now exp be a regular exponential map. As above, we split 
the t-integral in the second line below into two parts to get 

Nf(p) = J K*(p,0)Xf(p,O)d* p (6) 

(5.8) =/ [ ^( P: 8)K(exp p {t6),ex Pp (t8)) f(ex Pp (t6))dtd<j p (6) 

Js p M J 

= f W(p,v)f(ex Pp (v))dVo\(v), 
Jt p m 

where 

W = \v\- n+1 (j(p,v/\v\)K(ex Pp (v),exp p (v)/\v\) 

(5.9) V 

+ K*{p, -w/|«|)K(exp p (v), -e±p p (v)/\v\)j. 

Theorem 5.2. Let cxp p (w) satisfy (Rl) and (R4) and assume that for any (p,9) 6 suppK^, t6 is not a 
conjugate vector at p for t such that exp p (t9) € supp/. Then N is a classical ^>DO of order —1 with 
principal symbol 

(5.10) a p (N)(x,0 = 2w [ 5(Z(e))(K*K)(x,6)da x (6), 

JS X M 

where £(0) = (^O 3 , and S is the Dirac delta function. 

Proof. The theorem is essentially proved in Section 4 of [6], where the exponential map is related to a 
geodesic like family of curves. We will repeat the arguments there in this more general situation. 
Notice first that it is enough to study small enough Fix local coordinates x near p . By (R4), 

exp x {t&) — x + tm(t, 9; x), m(0, 9; x) = 9, 

with a smooth function m near (0, 9o,po). Introduce new variables (r,co) £ R x S X M by 

r = t\m(t,9;x)\, u — m(t, 9; x)/\m(t, 6;x)\, 

where | • | is the norm in the metric g{x). Then (r,oj) are polar coordinates for exp x (t9) — x = ruj with r 
that can be negative, as well, i.e., 

exp x (t9) — x + raj. 
The functions (r,ui) are clearly smooth got \t\ -C 1, and x close to po- Let 

J(i, 9;x) = det dt, v ( r , 

be the Jacobi determinant of the map (t,v) H> (r, cj). By (R4), J\t=o = 1, therefore that map is a local 
diffcomorphism from (— e,e) x S X M to its image for < £ <?C 1. It is not hard to see that for < e <€. 1 it 
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is also a global diffeomorphism, because it is clearly injective. Let t = t(x,r,ui), 9 — 6{x,r,oS) be the inverse 
functions defined by that map. Then 

t = r + 0(\r\), 9 = uj + 0(\r\), e*p(t0) = u + 0(\r\). 



Assume that the weight k in (2.2 1 vanishes for p outside some small neighborhood of po- Then after a change 
of variables, we get 

Nf(x)= / / A(x,r,ui)f(x + rw)drda x (ui), 
JS X M J 

where 

A(x, r, uj) = k^(x, 9(x,r, lu))k(x + ruj, uj + rO(l))J~ 1 (x, r, uj) 
with J as before, but written in the variables (x,r, uj). By [51 Lemma 2], AT is a classical ^DO with a 
principal symbol 



(5.11) 2tt / 5{£(u;))A(x,0,u))d(j x (u)) = 2tt / S^uj^k^x, uj)k{x, uj) da x {uj). 

Js x M Js„M 



□ 



Remark 5.3. Formulas (5.2) and (5.8) are valid regardless of possible conjugate points. In our setup, the 



supports of k, K* guarantee that exp Jt9), for (p, 9) close to (p , 9 ) reaches a conjugate point for t > but 



not for t < 0. Therefore, near the conjugate point q of p, the second term on the r.h.s. of (5.3), and (5.9), 
respectively, vanishes. 

6. The Schwartz kernel of N near the conjugate locus S 

We will introduce first three invariants. Let F : M — > N be a smooth orientation preserving map between 
two orientable Ricmannian manifolds (M,g) and (N,h). Then one defines detd-F invariantly by 

(6.1) F*(dVoljv) = (detdF)dVol M , 

see also [T^l X.3]. In local coordinates, 



. . det h(F(x)) , dF(x) 

We choose an orientation of S(po) near vq, as a surface in T Po M by choosing a unit normal field so that 
the derivative of detdexp po (t>) along it is positive on S(p). Then we extend this orientation to S(p) for p 
close to po by continuity. On Figure [2j the positive side is the one below S(p), if v is the first conjugate 
vector along the geodesic through (p, v). Then we choose an orientation of E(p) so that the positive side is 
that in the range of exp p . On Figure [2j the positive side is to the left of £(p). The so chosen orientations 



conform with the signs of £™ and y n in the normal form (3.2). 

Next we synchronize the orientations of T p M and M near q by postulating that exp p is an orientation 
preserving map from the positive side of S(p), as described above, to the positive side of £(p). 

For each p e M, the transformation laws in TT p M under coordinate changes on the base show that T p M 
has the natural structure of a Ricmannian manifold with the constant metric g{p). Then one can define 
det dexp p invariantly as above. Let d Vol p be the volume form in T p M , and let d Vol be the volume form in 
M. Then det d exp p is defined invariantly by 

(6.3) exp* d Vol = (det d exp p ) d Vol p . 
In local coordinates, 

a t-A /det 5 (exp v) 8 

detd6Xp ^ = \/ det g(p) det d- v ex ^ 
where, with some abuse of notation, g(p) is the metric g in fixed coordinates near a fixed p$, and g(exp p v) 
is the metric g in a possibly different system of fixed coordinates near go = ex P Po v °' 

(6.4) A(p,v) := |ddetdexp (v)\. 
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Since detdexp J9 (u) is a denning function for S(p), its differential is conormal to it. By the fold condition, 
A 0. One can check directly that A is invariantly defined on E. 

By (3.3), for (jp, v) £ S, the differential of exp p maps isomorphically T v S{p) (equipped with the metric on 
that plane induced by g(j>)) into T g E, with the induced metric. Let D be the determinant of exp p i.e., 

(6.5) D := det (dexp p \t v s{ p )) , 

defined invariantly by (6.f I. We synchronize the orientations of S(p) and S(p) so that D > 0. 

We express next the weight W(p,v) restricted to S in terms of the variables (p,q). For (p,q) £ S, 
v = exp~ (q), where we inverted exp p restricted to S. Let w = w(p, q) be defined as in (3.1 ) with v as above. 



Then we set, see also (15. 91), and Remark 5.3 



(6.6) W s (p,q) := W r ((p,ejt P - 1 (g))| E = |«| 1 - n /s ll (p ) t;/|«|)/s(g,-«;/|«|) 

For p close to po, E(p) divides M in a neighborhood of go m to two parts: one of them is in the range of 
exp p (w) for v near vo, that is the positive one w.r.t. the chosen orientation; the other is not. Let z'(p,q) be 
the distance from q to E(p) with a positive sign in the first region, and with a negative sign in the second 
one. Then for a fixed p, z' = z'(p, q) is a normal coordinate to £(p) depending smoothly on p, and S is given 
locally by z' — 0. Then z' is a defining function for E, i.e., S = {z' = 0} and d p . q z' ^ because d q z' ^ 0. 
Let z" — z"(p,q) £ R 2n_1 be such that its differential restricted to TT, is an isomorphism at (po,qo). Since 
dz" and dz' are linearly independent, z — z(z',z") are coordinates near (po,qo)- One way to construct z" 
is the following. Choose (z„ +1 , . . . , z 2n ), depending on p only, to be local coordinates for p, and to choose 
[z 1 ', 22, ... , z n ), depending on p and q, to be semi-geodesic coordinates of q near £(p). 

The next theorem shows that near S, the operator iV has a singular but integrable kernel with a conormal 
singularity of the type l/^/z*. 

Theorem 6.1. Near S(p), the Schwartz kernel N(p,q) of N (with respect to the volume measure) near 
(Po,1o) is of the form 

(6.7) N = Wv^L=(l + ^R(V?,z")), 

VADz' 

where Ws = Ws(z"), A = A{z"), D = D(z"), and R is a smooth function. 



Proof. We start with the representation (5.8). We will make the change of variables y = exp (v) for (p,v) 



close to (po,vo) as always. Then y will be on the positive side of S(p), and the exponential map is 2-to-f 



there. We split the integration in (5.8) in two parts: one, where v is on the positive side of S(p), that we 



call N + f, and the other one we denote by iV_/. Then 



(6.8) N±f(p)= Wf(y)(detdexp±(v)) dVo%), 

JS„M ■ 



where W is as in (6.6) but not restricted to S, and (exp^ ) _1 there is the corresponding inverse in each of 
the two cases. 

To prove the theorem, we need to analyze the singularity of the Jacobian determinant det dexp p (u) near 
S(p). ft is enough to do this at (po^o)- 

Let y — (y',y n ) be semi-geodesic coordinates near E(go), <7o = exPpo^o); and let yo correspond to qo. We 
assume that y n > on the positive side of S(p). In other words, y n = z'(po,q)- 

We have 

d Vol(y) = det (d„ exp p (w)) d Vol(w) 

The form on the left can be written as dVolz( p )(y') dy n ; while the one on the right, restricted to S(p), equals 
d Volg(p) (i/) dv n in boundary normal coordinates to S(p), where v n > gives the positive side of S(p). On 
the other hand, by (6.5), 

dVol s(p) (y')=DdVol s(p) (v'). 

We therefore get 

Ddy n =det (dexp p (v)) dv n . 
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By the definition of A, we have 
(6.9) 

Therefore, 

Since y n = for v n = 0, we get 



detd„exp («) = Av n (l + 0(v n )). 



Ddy n = A{l + 0(v n ))v n dv r 



y" 



(v 



n\2 



A 
2D 



(l + OK)). 



Solve this for v n and plug into ( |6.9[ ) to get 
(6.10) 



detdexp p (w) = ±^/2ADy n (l + 0±( V / F r )) • 



Here 0±( 
at « n = 



(6.8) 



y"j denotes a smooth function of ^/y™ near the origin with coefficients smooth in y' , that vanishes 
The positive/negative sign corresponds to v belonging to the positive/negative side of S(p). By 



(6.11) 



N±f(p) = / Wf(y) 



(l + ± (Vy n ))dVo\(y). 



y/2ADy r ' 

We replace A$, Dq in (6.11) by their values at y n = 0; the error will then just replace the remainder term 
above by another one of the same type. Similarly, W — W(p, v), where exp p (u) = q. Solving the latter for 
v = v(p, q) provides a function having a finite Taylor expansion in powers of s/y™ of any order, with smooth 
coefficients. The leading term is what we denoted by Ws that is a smooth function on E. 



With the aid of (6.2 ), it is easy to see that (6.11 ) is a coordinate representation of the formula (16. 71 at the 



so fixed p. When p varies near pq, it is enough to notice that since we already wrote the integral in invariant 
form, y n then becomes the function z'(p, q) introduced above. For z" we then have z"(p, q) = (x(p), y'(p, <?)). 
Finally, we note that another choice of z" so that (z' , z") are coordinates would preserve (6.7) with a possibly 
different R. □ 

7. N as a Fourier Integral Operator. Proof of Theorem 12.11 



We are ready to finish the proof of Theorem 2.1 By Theorem 6.1 near £, the Schwartz kernel of N has a 



conormal singularity at £, supported on one side of it, that admits a singular expansion in powers of W z[ 



+ ■ 



z' + . The Fourier transform of the latter is 



,-Lt/4(£-1/2 



+ iC 1/2 ) 



with a leading singularity 1/ A 
(7.1) 

where C+ = max(C,0), £- = (—£)+. The singularity near £ = can be cut off, and we then get a symbol of 
order —1/2, depending smoothly on the other 2n — 1 variables. Therefore, near E, the kernel of N belongs 
to the conformal class 7~"/ 2 (Af x M,S;C), see e.g., 18.2]. It is elliptic when K B (p , o )K i {q o , -w ) ^ 



by (5.9), (6.6). Therefore, the kernel of N near £ is a kernel of an FIO associated to the Lagrangian T*E. 
Moreover, the amplitude of the conormal singularity at £ is in the class S^^ 2 ' 1 ^ 2 (polyhomogeneous of order 
— 1/2, having an asymptotic expansion in integer powers of ICI 1 ^ 2 ); see a ls° (9.13) and (9.14). 



8. The two dimensional case 

Theorem 8.1. LetdimM — 2. Assume that (Rl) - (R5) are fulfilled. Then A/"*£ \ 0, near (po, £o, <7o, Vo)> 
is the graph of a local diffeomorphism T* M \ € (jp, £) <-> (q, rf) € T*M\ 0, homogeneous of order one in its 
second variable (a canonical graph). 

Proof. For (p, £) near (po,S,o), there are exactly two smooth maps that map £ to a unit normal vector. We 
choose the one that maps £o to vo/l^ol- Then we map the latter to v € S(p). Since the radial ray through 
v is transversal to S(p), that map is smooth. Knowing v, then we can express q — exp p (v) g £(p) and 
-exp p (u) as smooth functions of (p,£) as well. Then in local coordinates, r\ = £jdexp* (w)/dq, see 



w 



(4.12), that in particular proves the homogeneity. 



By (R5), this map is invertible. 



□ 
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Figure 3. The 2D case 



The principal symbol of X*X in the geodesies case, see Theorem 5.1 and (5.5), is given by 
(8.1) a p (X*X)(x,0 = ^\<x,^/\^\)\\ 

where £- L is a continuous choice of a vector field normal to £ and of the same length so that at p — po, 
£o~/l£o~l = 0o, — £o"/l — £o~l = 0o! therefore, the sign of the angle of rotation is different near £ and near — £ - 



Notice that (5.5) in the two dimensional case is a sum of two terms but we assumed that k is supported neat 



(poj0o)j therefore only one of the terms is non-trivial. A similar remark applies to (5.10). 



Theorem |6.1| takes the following form in two dimensions, in the Riemannian case. 



Corollary 8.1 

of Theorem \6.1 

(8.2) 
where 



Let n = 2 and let exp be the exponential map of a Riemannian metric, 
we then have 

V2 



With the notation 



N = W s ^^(l + Vz'R(Vz',z")) 



D 



is evaluated at v € S(p) such that q = exp (v), 



-r^det d exp p (v) 

and d/dN stands for the derivative in the direction of N p (v). 



Proof. Note first that B ^ by the fold condition. Let <fi be the (acute) angle between S(p) and N p (v) at 
v. Since N p (v) is orthogonal to the radial ray at v, we can introduce an orthonormal coordinate system at 
v with the first coordinate vector being i>/|v|, and the second one: the positively oriented unit vector along 
N p (v), that we call £. Let us parallel transport this frame along the geodesic 7 P ,„; and invert the direction of 
the tangent vector to conform with our choice of w at q. In particular, this introduces a similar coordinate 
system near the corresponding vector w at q in the conjugate locus. In these coordinates then 

3.3) dexp p (u) 







3/\v\ 



where j is uniquely determined by J(t) = j(t)E(t), where J(t) is the Jacobi field with J(0) = 0, J'(0) = £, 
and H(t) is the parallel transport of £, compare that with ( |4.1| . The extra factor l/\v\ comes from the fact 
that we normalize v now in our basis, so that the result would be the Jacobian determinant. Then the Jacobi 
determinant detdexp p (u) is given by —j/\v\. In particular, for (p,v) £ S we have dexp p (u) = diag(— 1,0). 
Note that j depends on v as well, therefore its differential that essentially gives ddetdexp p (zj) depends on 
the properties of the Jacobi field under a variation of the geodesic. 
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Now, it easily follows from the definition (6.5| of D that 

D = sin0. 

On the other hand, ddetdexp p (w) is conormal to S(p), therefore, the derivative of detdexp p (u) in the 
direction of N p (v) satisfies 



-jjy det d exp p (u) = |ddet dexp p (w)| sine/) = Asmtj) = AD. 



□ 



9. Resolving the singularities in the geodesic case 

Let, as before, (po;<7o) be a pair of fold conjugate points along 70, and X be the ray transform with a 
weight that localizes near j . We want to see whether we can resolve the singularities of / near p and near 
qo knowing that X / G C°°, and more generally, whether we can invert X microlocally. Assume for simplicity 
that p ^ q a . 

We will restrict ourselves to the geodesic case only but the same analysis holds without changes to the 
case of magnetic geodesies as well. We avoid the formal introduction of magnetic geodesies for simplicity of 
the exposition. Assume also that 



(9-1) 



n(p, 9)n{q, -w/\w\) ^ 0, for (p, 6) € U , 



where (<?, w) are given by ( |3.1[ ), and Wg)Wo3 (po, 6>o). This guarantees the microlocal ellipticity of the ^PDO 
A near Af*(po,vo) and N*(qo, w ) in Theorem |2.1| see Theorem 5.1 



9.1. Sketch of the results. We explain the results before first in an informal way. As we pointed out 
in the Introduction, Xf(j) for geodesies near 70 can only provide information for WF(/) near J\f* 70, and 
does not "see" the other singularities. The analysis below based on Theorem |2.1| shows that on a principal 
symbol level, the operator \D\ X / 2 F behaves as a Radon type of transform on the curves (when n = 2) or the 
surfaces (when n > 3) £(p). Similarly, its adjoint behaves as a Radon transform on the curves/surfaces £(g). 
Therefore, there are two geometric objects that can detect singularities at p conormal to v : the geodesic 
7o = 7p ,uo (and those close to it) and the conjugate locus £(<7o) through po (and those corresponding to 
perturbations of v ). We refer to Figure |4j 

When n = 2, the information coming from integrals along the two curves (and their neighborhoods) may 
in principle cancel; and we show in Theorem |9.2| that this actually happens, at least to order one. When 
n > 3, the Radon transform over T,(q) 9 p competes with the geodesic transform over geodesies through 
p. Depending on the properties of that Radon transform, the information that we get for ±£0 may or may 
not cancel because £0 is conormal both to 70 and £(<7o)- On the other hand, for any other £1 conormal to 
vo but not parallel to £o> the geodesic 70 (and those close to it) can detect whether it is in WF(/) but the 
Radon transform restricted to small perturbations of v (and therefore of go) wm n °t- Thus, we can invert 
N microlocally at such (pojCi)- 

Now, when n > 3, we may try to invert N even at £0 by choosing u's close to vq but normal to £o- If 
£0 happens not to be conormal to the corresponding conjugate locus T,(q(p ,v)) at po, we can just use the 
argument above with the new v. In particular, if the map (4.3) is a local diffeomorphism, this can be done. 



This suggests the following sufficient condition for inverting N at (po,Ci)' 



(9.2) 



3#i G S Po M, so that /c(po>#i) 7^ 0, £i(#i) = 0, and £1 is not conormal to T,(q(p ,9i)) at Pq. 



Above, £(q(po,0i)) is the conjugate locus to the point q that is conjugate to po along 7p ,9 1 . We normally 
denote that point by q(po,vi), where V\ G S(po) has the same direction as 0\. 



In case of the geodesic transform, one could formulate (9.2 1 in terms of the map (4.3) as follows: 
(9-3) 



3vi G S(po), so that k(j>q, v±/\vi\) 7^ 0, £1(^1) = 0, and £1 is not the image of v\ under the map (4.3) at po 



Section 10.4 we present another example, where (9.2) fails 



In Section 10.3 we present an example where (4.3) is a local diffeomorphism, therefore (9.2) holds. In 
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9.2. Recovery of singularities in all dimensions. We proceed next with analysis of the recovery of 
singularities. 

Let xi,2 be smooth functions on M that localize near and qo, respectively, i.e., supp Xi C U\, supp X2 C 
U 2 , where U± l2 are small enough neighborhoods of p and q, respectively. Assume that xi, Xi equal 1 in smaller 
neighborhoods of po, qo-, where /i, f 2 are supported. Then / :— f\ + f 2 is supported in Ui U U2 and we can 
write 

(9.4) XiNf = A 1 f 1 + F 12 f 2 , 

where Ai — Xi^Xi is a ^DO by Theorem |5.2| while F\ 2 = Xi^X2 is the FIO that we denoted by F in 



Theorem 2.1 By (R5), we can do the same thing near qo to get 
(9.5) X2Nf = A 2 f 2 + F 21 f u 

where A 2 = X2NX2, F 2 i — X2^xi- It follows immediately that F 2 i — F* 2 . Recall that F12 — F in the 
notation of Theorem |2.1| Assuming X*Xf 6 C°°, we get 



(9.6) Ax h + Ff 2 e C°°, A 2 f 2 + F*h e C°° . 
Solve the first equation for / 2 , plug into the second one to get 

(9.7) (Id - A 2 x F*A^F)f 2 e C°° near (q , ± Vo ) , 

where A^ 1 , A^ 1 , denote parametrices of A\, A2 near (po,±£o)> an( i (90;i?7o)j respectively. The operator 
in the parentheses is a ^PDO of order if the canonical relation is a graph, that is true in particular when 



2, by Theorem 8.1 In that case, if Id — A 2 ~ 1 F*A 1 ~ 1 F is an elliptic (as a ^PDO of order 0), near (q , ±%), 
then we can recover the singularities. Without the canonical graph assumption, if it is hypoelliptic, then we 
still can. 

Another way to express the arguments above is the following. Since xi,2 together with k restrict to conic 
neighborhoods of (po ± £0), and (qo ± 770), respectively, and ^1,2, F, F* have canonical relations of graph 
type that preserve the union of those neighborhoods, we may think of / = /1 + / 2 as a vector / = (/ 1; / 2 ), 
and then 

(9-8) F=(p I 

The operator Id — A 2 ~ 1 F*A 1 ~ 1 F can be considered then as the "determinant" of F, up to elliptic factors. 
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Theorem 9.1. Let the canonical relation of F be a canonical graph. With the assumptions and the notation 
above, if the zeroth order ^ DO 

(9.9) Id-A^ 1 F*Ai 1 F 

is elliptic in a conic neighborhood of (qo,±r]o), then Xf 6 C°° near (pq,9q) (or more generally, Nf € C°° 
near po and qo ) implies f <E C°° . 

In the geodesic case in two dimensions, the principal symbol of F*A^ F is always 1, see the Propo- 
sition 19.11 below. 



When n > 3 and F is of graph type, then A 2 1 F*A 1 1 F is of negative order, therefore we can resolve the 
singularities. 

Corollary 9.1. Let n > 3 and assume that the canonical relation of F is a canonical graph. Then the 



conclusions of Theorem 9.1 hold, i.e., Xf € C°° near (po,0 ) (or more generally, Nf € C°° near p , q ) 
implies f <E C°°. 

Proof. In this case, A^F is an FIO of order 1 — n/2 with the same canonical relation is F. Similarly A 2 1 F* 
is an FIO of order 1 — n/2 with a canonical relation that is a graph of the inverse canonical map. Their 
composition is therefore a ^DO of order 2 — n < 0. Its principal symbol as a ^DO of order is zero. The 
corollary now follows from Theorem |9.1| □ 

In Section |10.3| we give an example where the assumptions of the corollary hold. Note that those 
assumptions are stable under small perturbations of the dynamical system. 



When the graph condition does not hold, the analysis is harder. Then (4.3) is not a local diffeomorphism 



If its range is a lower dimension submanifold, for example, we can at least recover the conormal singularities 



to 6*o away from it, as the corollary below implies. Note that below, (b) implies (a). Also, (9.1 ) is not needed; 
only ellipticity of k at (poj^o) suffices. 

Corollary 9.2. Let Xf E C°° for 7 near 7o . Then 

( a ) If £1 € Tp M \ is conormal to Vq but not conormal to £(<fo) (not parallel to £,0), then 

(pd,6)0WF(/). 



(b) The same conclusion holds if condition (9.2) or the equivalent (9.3) is fulfilled. 



Proof. Note first that A\ is elliptic at (j»qjC) by (9.1) and Theorem |5.1[ b). By the first relation in (9.6), 



(fOiCl) € WF(/i) if and only if (p<j,£i) E WF(F/ 2 ). To analyze the latter, we will use the relation 
WF(F/ 2 ) C WF'(F) o WF(/ 2 ), see % Thm 8.5.5]. Note also that in the notation in Thm 8 5.5] , 



W¥(F)x is empty. By Theorem 6.1 WF'(F) consists of those points in the canonical relation C, see (2.5), 
for which the conormal singularity in ( |6.7[ ) is not canceled by a zero weight. 

Now, let £1 be as in (a). Since £1 is separated by ±£ by a conic neighborhood, one can choose a weight 
X on SM that is constant along the geodesic flow, non-zero at (po,9q) and supported in a flow-out of a 
neighborhood V of it small enough such that the conormals to the corresponding conjugate loci at po stay 



away from a neighborhood of £1. In the geodesies case, the condition is that the map (4.3) restricted to V, 
does not intersect a chosen small enough conic neighborhood of ±£o- This can always be done by continuity 
arguments. Then left projection of WF'(F) will not be singular at (po,£i), and therefore, Ff 2 will have the 
same property regardless of the singularities of / 2 . 

Statement (b) follows from (a) by varying v near vq in directions normal to £1. □ 



9.3. Calculating the principal symbol of (9.9) in case of Riemannian surfaces. Let exp be the 

exponential map of g, and let n > 2. We will take n — 2 late r. Recall that the leading singularity of 
the kernel of N near E is of the type (z^)^ 1 / 2 , by Theorem 6.1 We will compose F with a certain 'J'DO 



R so that this singularity becomes of the type S(z'). Then modulo lower order terms, FRf{p) will be a 
weighted Radon transform over the surface In 2D, that will be an X-ray type of transform. We are 

only interested in this composition acting on distributions with wave front sets in a small conic neighborhood 
W of {q ,±r]o). 
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The Fourier transform of (z' + ) 1 / 2 is given by (7.1). Its reciprocal is 



( MC)C l/ 2 _ iM _ C)( _ C) l/ a ) = 7r -l/ 2e W4( MC) _ i/ l (-C))|C| 1/2 , 

where h is the Heaviside function, and \(\ is the norm in T*M . We fix p near p and local coordinates 
x = x(p) there, and we work in semi-geodesic coordinates y — y(p, q) near qq normal to £(p) oriented as in 
section [6j Let x denote local coordinates near q . Let R be a properly supported ^DO of order 1/2 with 
principal symbol, equal to 

(9.10) r( 2 /,7 7 ) = 7r- 1 / 2 e-/ 4 (% n )-^(-^))|r?r /2 ro(y ) 77) I 

in W, outside some neighborhood of the zero section, where r$ is a homogeneous symbol of order 0, an even 
function of rj. Note that 

(9.11) \ r \ 2 =^\ v \rl 

The appearance of the Heaviside function here can be explained by the fact that iV*E has two connected 
components: near (pa, qo, —£o, 770) and near (po, 9o, ?o> — ^0); and the constants needs to be chosen differently 
in each component. 

We start with computing the composition 

(9.12) FR. 



Since the kernel of (9.12) is the transpose of that of RF 1 , we will compute the latter; and we only 
need those singularities that belong to W. Denote by F(p, q) the Schwartz kernel of F. Then the kernel 
F'(q,p) = F(p,q) of F' (with the notation convention F'f(q) — J F'(q,p)f(p) d Vol(p)) can be written as 
F'(q(x, y),p(x)) that with some abuse of notation we denote again by F'(y,x). Then 

(9.13) F'(y, x) := (27T)" 1 J ^F'(y', Vn ,x) d Vn , 

where F' is the partial Fourier transform of F w.r.t. y n , and there is no summation in y n r\ n . By Theorem 



6.1 



and fl7.1| ), 

(9.14) F'(y', Vn , x) = ■K 1 ' 2 e-™' A (h{i ln ) + ih(- Vn )) |r 7n |- 1 /2 G ( S) y ', r, n ) 

where G is a symbol w.r.t. rj n , smoothly depending on (x,y') with principal part 

V2 



G :=W S - 



'AD 



Moreover, by Theorem 6.1 G has an expansion it terms of positive powers of \r) n \ 1 ^ 2 . In particular, G — Go 
is an amplitude of order —1/2 that contributes a conor mal distribution in the class I^ n / 2 ^ 1 / 2 (M x M, E; C), 
see, e.g., [H Thm 18.2.8]. By the calculus of conormal singularities, see e.g., [U Theorem 18.2.12], the 
kernel of FR is of conormal type at y n = as well, with a principal symbol given by that of F multiplied 
by r\ y n = o ir] >=o- That principal symbol coincides with the full one modulo conormal kernels of order 1 less 
that the former, see the expansions in [5] preceding Theorem 18.2.12. Since we assumed that ro is an even 
homogeneous function of 77 of order 0, ro(y' , 0, 0, rj n ) is a function of y' only for rj in a conic neighborhood of 
(0, ±1), equal to r(y, 0, 0, 1). Therefore, the principal part of r(y, D y )F'(-, x) is 

(9.15) (27T)- 1 / e i ^"G ( a; ,77')ro(y',0,0,l)d7 7n = W^^=r (y', 0,0, l)S(y n ), 

J VAD 

and the latter is in I^ n / 2+1 / 2 (M x M, S;C). The "error" is determined by the next term of the principal 
symbol of the composition FR with G replaced by Go, that is of order 1 lower and by the contribution of 
G = Go that is of order —1/2 lower. Since the coordinates (y',y n ) depend on p, as well, ro(y', 0, 0, 1) is 
actually the restriction of ro to J\f*Yl(p). So we proved the following. 



THE GEODESIC X-RAY TRANSFORM WITH FOLD CAUSTICS 21 



Lemma 9.1. Let r be as in Then modulo I- n l 2 {M x M, E; C), FR G jV2-n/2(M x M,£;C) 

reduces to the Radon transform 

FRf{p) ~ / a/dS, a := ^ S ^L, 

where dS is the Riemannian surface measure on £/ia£ we previously denoted by dVoljv p ). 

In two dimensions, this is an X-ray type of transform. In higher dimensions, this is a Radon type of 
transform on the family of codimension one surfaces S(p). 
In what follows, n = 2. 
We will compute RF*FR next. We have 

(9.16) [ FRfFRhdVol- [ [ (af)(z')dS{z') [ (ah)(q) dS{q) d Vol(p) 

modulo terms of the kind (P/, /i), where P is a ^PDO of order —3/2 or less. 

In the latter integral, p parameterizes the curve S(p), while q G S(p) parameterizes a point on it. Another 
parameterization is by p and £ G S*M with £ oriented positively; then q — exp p (w), where v G and 
= 0. For the Jacobian of that change we have 

(9.17) dS(q)dVo\(p) = DdVo\ s(p) (v)dVo\(p) = ^ da p (0 d Vol(p), 

and we recall that dcr p denotes the surface measure on S p M, that in this case is a circle. The canonical map 
(p,£) — > (9,77) is symplectic, and therefore preserves the volume form dpd£. Set 

(9.18) K:=\ V (p,0\m- 

Then this map takes S*M into {(q,r)) G T*M; \rj\ — K}. Project that bindle to the unit circle one, and set 
fj = T]/\t)\. Then we have the map (p, £) —> (q, fj), and d Vol(p) da p (£,) — K 2 d Vol(g) da q {i)). 
When we perform those changes of variables in (9.16), we will have 

(9.19) dS(q) d Vol(p) = vol(g) da q {n), 

COS 

where p G M, q G (q,r]) G S*M, and we removed the hat over 77. Let w is the corresponding 

vector in S(q) normal to 77. That parameterizes the curves over which we integrate by initial points 
q and unit conormal vectors 77. The latter can be replaced by unit tangent vectors w = w/\w\; then 
d Vol(g) dcr q (ri) — d Vol(q) da q (w). Let us denote the so parameterized curves by c q ^(s), where s is an 
arc-length parameter. 

It remains to notice that the integral w.r.t. z' G is an integral w.r.t. the arc-length measure on 
that we denote by s. Then performing the change of the variables (p, q, z') 1— >• {q, w, z') in ( 9.16 ), we get 

(9.20) [ FRfFRhdVol- [ (af)(c q ^(s))a(q,-w)h{q)ds do - g (^) d Vol(g). 

J JRxS„MxM ' COS0 



Therefore, we get as in (5.2 1, (5.4 1 

y/dct(g(q)) J cos0 p(q,g') 

K -=4= t ]nW .^ lwd ^ ML dmW) 

^det(g(q)) J 1 y n Acoscj) p(q, q') 



For the directional derivatives of detdexp p (7j) = —J'/\v\, see (8.3), we have that the derivative along the 
radial ray is |J'(1)|/|«| by absolute value, while the derivative in the direction of S(p) vanishes. That implies 

Acos(t>= |J'(1)|/M =K/\w\. 
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Therefore, 
(9.22) 



R* F* FRf{q) 



Vdet(.g(g)) 



2K\r \^ {p) \ 2 \Wx\ 2 \w\ 2 ^^dVo\( q >). 



Here (p, v) is defined as follows. It is the point in SM that lies on the continuation of the geodesic through q, 
q' to its conjugate point near p$, The weight k restricts q' to a small neighborhood of 70. Next, A2 restricts 
q' near qg. 

We compare (9.221 with (5.4) and (5.5). Notice that the Jacobian term in (5.4) at the diagonal equals 
-y/det g and therefore cancels the factor in front of the integral in the calculation of the principal symbol. We 
therefore proved the following. 

Lemma 9.2. Let n = 2. Then R*F*FR is a ^DO of order —1 with principal symbol modulo 5~ 3 / 2 at (q, rf) 
near {qo,r/o) given by 

4^|7 ? |- 1 |roU. s(p) | 2 |K(p,«/|i;|)| 2 |K(g,-«;/| W |)| 2 

Here w/\w\ is a continuous choice of a unit vector normal to r/ at q, so that (q,w/\w\) = (qo,wo/\wo\) when 
(QtV) = (qotVo): an d v /\ v \ * s a parallel transport of — w/\w\ from q to its conjugate point p along the geodesic 

Later we use the notation w = n- L /\r]' L \, and v = £ /|£ |. 



Proposition 9.1. Let n = 2. Then 
is a ^DO of order -1/2. 



Id— A 2 1 F*A^ 1 F 



1-1/2 



Proof. We apply Lemma 9.2 with 7r 1//2 e 17r / 4 |?7| 1 / 2 ro being the principal symbol of A 2 1//2 , see (9.10), where 



is a parametrix of A^* near (q , ±f?o)- To this end, choose 



see (8.1 ). Note that n{q, w/\w\) — k(p, —v/\v\) = because of the assumption on supp k. Then |ro|AA*s( P ) 

2~ 1 l 2 \K(q, — w/|it;|)| _1 , where w is as in (3.1 ). The choice of ro yields RR* = A 2 mod So Lemma 

implies that R*F*FR, and therefore RR*F*F and A 2 ~ 1 F*F, have principal symbol 



9.2 



a p (A 2 1 F*F)(q, V ) = 2ttK| K 0,£ x /I£- 



We only need to insert A 1 1 between F* and F. By (TUJ Thm 25.3.5], modulo ^DOs of order 1 lower, 
the principal symbol of A 2 ~ 1 F*A^ 1 F is given by that of A 2 ~ 1 F*F multiplied by the principal symbol 
(2ir\n(p, w)| 2 |/|£|) 1 °f pushed forward by the canonical map of F. In other words, 



a p (A 2 - 1 F*A^ 1 F)(q,r 1 ) 



2n\K(p,e/\^\ 



'-K[2n\K((p,e/\z x \)\ 2 /m,n)\)] 1 



□ 



The following lemma is needed below for the proof of Theorem 9.2 



Lemma 9.3. Let n\ and k both satisfy the assumptions for k in the Introduction, and let n(po, #0) 7^ 0. Let 
XG $° have essential support near (po, ±£0) U (qo, ±?7o) and Schwartz kernel in (Ui x U±) U (U2 x Uz)- Then 
there exists a zero order classical ^ DO Q with the same support properties so that 

QX*X kX = X* Ki X kX , modr^ 2 (M x M, AuA^*S, C), 

where A is the diagonal. In particular, QX*X K x — X*X K x '■ H s — > H s + 3 / 2 i s bounded for any s. 
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Proof. We define Q = Qi + Q2 where Qi,2 have Schwartz kernels in U\ x U\ and XJ% x [7 2 , respectively. 
Following the notation convention in (9.8), Q = diag(Qi, Q 2 )- 
Then we choose Q\ to have principal symbol 

(9.23) Kx(jp,^/\e\)/-<p,e\/\^) 

in a conic neighborhood of (pci)±£o) with the same choice of £~ L as in 
principal symbol 



.11. Next, we choose Q2 with a 



(9.24) 

in a conic neighborhood of (qo,±r]o) 



^iiq,v ± /\v x \)/^{Q,v ± \ 

Then 



QX* K x R = 



Q1A1 
Q2F* 



QiF 

Q2A2 



Then, see (8.1), 



o-p{QlM) = 27r(R lK )(p^ ± /\^\)\, <j p (Q 2 A 2 ) = 2TT(R lK )(q,r 1 ± /\ V ± \). 

For QiF, Q2F* , we use the arguments used in the proof of Lemma [9. 1| A representation of the Schwartz 
kernel of F' as a conormal distribution is given by (9.13 1 . The composition Q2F* then is of the same conormal 
type with a principal symbol equal to the co mpl ex conjugate of th at o f F' multiplied by the symbol (9.24) 
restricted to 7V*£. This replaces fe" = R in (6.6) by Ri. Since in (6.6), «' = R we get that Q2F* is of the 



same conormal type with leading singularity as in Theorem |6.1[ with 

VF S = \v\~ 1 R(p,v/\v\)K 1 {q,-w/\w\). 

This is however the leading singularity of X2X* 1 X K \i- 

The proof for Q\F is the same with the roles of p and q replaced. 



□ 



9.4. Cancellation of singularities on Riemannian surfaces. Assume in all dimensions that there are 
no conjugate points on the geodesies in M, and that dM is strictly convex. Let Mi D M be an extension of 
M so that the interior of Mi contains M be as in Remark |5.2| Then if k 7^ 0, 

(9-25) \\f\\L*(M) < C\\X*Xf\\ H i {Ml) + C k \\f\\ H - HM) , V/gL 2 (M), 

for all k > 0, see [HI [6], and [17] for a class of manifolds with conjugate points. When we know that X is 
injective, for example when the weight is constant; then we can remove the H~ k term. The same arguments 
there show that for any s > 0, 

(9.26) WfWHHM) < C\\X*Xf\\ H s +HMl) + C k \\f\\ H - HM) , v/ g //,:;. \/K 

Consider Xf parameterized by points in d-\-SM±, that defines Sobolev spaces for Xf as in section 5.1 Then 

(9.27) ||/|| ff .(jio < C\\Xf\\ H . +1 , HB+SMl) + C k \\f\\ H - k(M) , V/ G H°(M), s > 0. 

Indeed, in Proposition |5.2[ one ca n com plete Mi a nd H to closed m anifol ds, and then we would get that X_ 
H s+1 / 2 is bounded. Then (9.27) follows by (9.26). Estimate (9.27) is sharp in view of Proposition 



5.2 



In the following theorem, we show that ( |9.25| , ( |9.27[ ) fail in the 2D case, with a loss at least of one derivative 
in the first one, and 1/2 derivative in the second one. 

Theorem 9.2. Let n = 2, and let 70 be a geodesic of g with conjugate points satisfying the assumptions in 
section^ Then for each f 2 G H S (M), s > 0, with WF(/ 2 ) in a small neighborhood of (q$,±rf ), there exists 
fi G H S (M) with WF(/i) in a some neighborhood o/(po,±£ ) so that 

Xf G H s+3/4 and X* Xf G H s+3/2 , where f := fi + fa. 

In particular, if (M, g) is a non-trapping Riemannian surface with boundary with fold type of conjugate points 
on some geodesies, none of the inequalities (9.25), (9.21) can hold. 

Remark 9.1. It is an open problem whether we can replace ff s + 3 / 2 and H s+2 above with C°°. See 
Section 10.1 for an example where this can be done. 
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Remark 9.2. If there are no conjugate points, one has Xf G £P +1/2 , X*Xf £ H s+1 . Therefore, the 
conjugate points are responsible for an 1/4 derivative smoothing for Xf, and an 1/2 derivative smoothing 
for X*Xf 

Proof. Let f 2 be as in the theorem. Set 

fx = -A?Ff 3 , 

where, as before, Aj 1 , A^ 1 are parametrices of Ai t2 in conic neighborhoods of (po,±£o) and (qo,±r]o), 
respectively. Then /i belongs to H s and has a wave front set in small neighborhood of (po ± £o)j by 
Theorem 2.1 By construction and by (9.4), 

(9.28) 



Next, by (9.28), 



XiX*Xf G C c 



A 2 f 2 + F*h = A 2 f 2 - F*A^Ff 2 = (A 2 - F*A^F)f 2 



The operator in the parentheses is a 'I'DO of order —3/2 by Proposition 9.1 Therefore, see (9.5), 

X2 X*Xf = A 2 f 2 + F*fi G H s + 3 ^ 2 . 

We therefore get X*Xf G £P +3 / 2 ([/i U t/ 2 ). 

To prove X/ G H s+1 , note first that above we actually proved that 



(9.29) 



X*X(ld - A- 1 F) X : H S {U 2 ) — -> H S+3 ' 2 {U 1 U C/ 2 ) 



is bounded, being a 'i'DO of order —3/2, where % denotes a zero order VPDO with essential support in a 
small neighborhood of (po>±f?o) and Schwartz kernel supported in U 2 x £/ 2 . 
Our goal is to show that 



X(Id-A^F) X :H s (U 2 



h3/4 



(H) 



is bounded. It is enough to prove that 



(9.30) 

for any *DO F 



X*(Id - A 1 - 1 P)*X*P 2s+3/2 X(Id - Ar X F)x : i/ s ((7 2 ) — > i?- s (C/ 2 ) 



(9.29) 



2s+3/2 



of order 2s + 3/2 on H. All adjoints here are in the corresponding L 2 spaces. By 



(9.31) Q 2s+3 /2***(Id - A- 1 F) X : H S (U 2 ) 

is bounded for any ^DO Q 2s +3/ 2 of order 2s + 3/2. 



H- S {U 2 ) 



To deduce (9.30) from (9.31 1, it is enough to "commute" X* with P 2s +3/ 2 in (9.30). Let 2s + 3/2 be a non- 



5.2 



we use the fact that X 



P 



2s+3/2 



2s+3/2 



X)*, and 



negative integer first. As in the proof of Proposition 

P 2s+2 X f is a finite sum of X-ray transforms with various weights of derivatives of / of order not exceeding 
2s + 2. Thus we can write 

(9.32) 



X*P- 



2s+2 



where Qj are differential operators on % of degree 2s + 3/2 or les s, a nd Xj are like X in (2.1) but with 
different weights still supported where k is supported. By Lemma 



9.3 



} 3 X*X 



RjX*X, where R 3 is a 



'I'DO of the same order as Qj. The proof of (9.30 1 is then completed by the observation that x*(Id — A 1 1 F)* 
maps continuously H s into itself, since the canonical relation of F is canonical graph. □ 

10. Examples 



In this section, we present a few examples. We start in Section |10.1| with the fixed radius circular 
transform in the plane, where we can have cancellation of singularities similarly to Theorem 9.2 but we 



show that this happens to any order. Then we consider in Section 10.2 the geodesic X-ray transform on the 
sphere, where the conjugacy is not of fold type, but a similar result holds. Next, in Section [10. 3[ we study 
an example of magnetic geodesies in the Euclidean space R 3 with a constant magnetic field. We show that 
then the canonical relation of F a canonical graph, and therefore, one can resolve the singularities. Finally, 



in Section 10.4 we present an example of a Riemannian manifold of product type where the graph condition 
is violated. 
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10.1. The fixed radius circular transform in the plane. Let R be the integral transform in R 2 of 
integrating functions over circles of radius 1. We fix the negative orientation on those circles; then for each 
(cc, £) € SR. 2 , there is a unique unit circle passing through x in the direction of 9. It is very easy to see that 
the first conjugate point appears at "time" ir. The next one is at 2w, that equals the period of the curve. If 
one originally chooses / supported near, say (0, 0) and (2, 0); and chooses 70 to be the arc of the circle that 
is a small extension of {\xi — 1| 2 + x\ = 1, X2 > 0}, then we are in the situation studied above. On the 
other hand, if we do not impose any assumptions on supp/, we will get contributions that are smoothing 
operators only. Therefore, we do not need to restrict supp/. 
The conjugate locus is the circle 

nx) = {y, \y-x\ = 2) 

that is the envelope of all circles of radius 1 passing through x. It follows immediately that 

S x (v) = {v; \v\ = tt}, N x (v) = Re ia (2/vr, -1), 

where we used complex identification to denote rotation by the angle a = arg(u). Hence, S is a fold conjugate 
locus. The other assumptions of the dynamical system are easy to check. 

It is much more natural to parametrize those circles by their centers, wc use the notation C(x). Then 



(10.1) Rf(x)= fd£= f(x + u)dt UJ = f(z + e ia )da. 

JC{x) J\u\=l JO 

Those circles are also magnetic geodesies w.r.t. the Euclidean metric and a constant non-zero Lorentzian force. 
Note that the "geodesies" are naturally parametrized by a point in R 2 as well (and that parametrization 
reflects the choice of the measure w.r.t. which we take R* , it is not hard to see that this is the same measure 
that we had before). 

10.1.1. R as a convolution. It is well known and easy to see that R is a convolution with the delta function 
<5gi of the unit circle 

Rf = S s i * f. 

Fourier transforming, we get 

(10.2) R = 2nF- 1 M\i\)F, 
where Jo is the Bessel function of order 0. This shows that 

(10.3) R*R = {2nfF- 1 Jl{\i\)F. 

Note that Jq (|£|) is not a symbol because it oscillates. In principle, one can use this representation to analyze 
R*R. 

10.1.2. Integral representation. We write 

(Rf, Rh)= I I f(x + ui) d£ u I h(x + 9) d£ dx 

(10.4) J J ^ 

f{y + co-9)h(y)de uj d£e dy. 



M=i J\e\=i 



Therefore, 

(10.5) R*Rf(x)= [ ( f(x + u + 9)d£ w d£ e , 

J\u\=i J\e\=i 



compare with (5.1) 
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We will make the change of variables z = oj + 9. For < \z\ < 2, there are exactly two ways z can be 
represented this way. Write uj = e la , 9 = e lf3 . Since dl^ — da, dig — d/3, and dzi A dz-i = (— 2z)~ 1 dz A dz, 
we get 

dzi Adz 2 = -^ (ie ia da + ie i/3 d/3) A (-ie~ i£ Ma - ie" i/3 d/?) = sin(/3 - a) da A d/3 
= sin(/3-a)d^Ad4. 

It is easy to see that \/3— a\ equals twice the angle between z = uj + 9 and 9. Let r = \z\. Then r/2 = cos ^ a ~^ . 
Elementary calculations then lead to 

T I 

sin \a — f3\ = - v 4 — r 2 . 



Therefore, (10.5) yields the following. 



Proposition 10.1. Let R be the circular transform defined above. Then 

(10.6) R*Rf(x)= [ —- L=f(y)dy, r:=\x-y\. 

Jr<2 rv4 — r £ 

10.1.3. R*R as an FIO. The kernel has singularities near the diagonal x = y, and also near 

S = {|x-y| =2}. 

That singularity is of the type (2 — \x — y|) -1 / 2 , and for a fixed x the expression 2 — \x — y\ measures the 
distance from the circle to the point y inside that circle. We therefore get the same singularity as in 
Theorem 16. II Note also that 

(10.7) Af*T, = {(x, x ± 2£/|£|, £, -?);^R 2 \ 0}. 

Based on Proposition |10.1| and Theorem 1 2. 1| we conclude that R*R is an FIO of order —1 with a canonical 
relation C of the following type. We have that (x, £, y,rj) E C if and only if (y, rj) = (x, £) (that gives us the 
fDO part), or (y, V ) = (x ± 2£/|£|,0- 

This can also be formulated also in the following form. 

Theorem 10.1. Let R be the circular transform defined above. Then, modulo , 

(10.8) R*R = A Q + F+ + F_, 

where Aq, F+ and _F_ are Fourier multipliers with the properties 

(a) A = 4tt|L>|- 1 mod v^ -1 ; 

(b) F± are elliptic FIOs of order —1 with canonical relations of a graph type given by 

(10.9) T ± :(x,d)y->(x±2Z/\t;lO- 

(c) F_=F* + . 

Proof. We start with the Fourier multiplier representation (10.2 1. The leading term of (2-7r) 2 Jq (|£|) is 

8tt 8tt / 2 e 2i l ? l e ~ 2i l€l 

(10.10) || cos 2 (|e - tt/4) = ^(1 + sin(2|e|) = 2tt + - 

Those three terms are the principal parts of the operators in (10.8|. The first one gives 47r|_D|~ 1 , while the 
second and the third one are FIOs with phase functions (j>± = (x — y) ■ £ ± 2|£|. A direct calculation show 
that the canonical relations of F± are given by (10.91, indeed. For the complete proof of the theorem, we 
need the full asymptotic expansion of Jo- 

We recall the well known expansion of Jq(z) for z —> oo: 



J (z) ~ y/2/(irz) (P(z) cos(z - tt/4) - Q(z) sin(z - tt/4)) , 

where 

oo oo 

P(z) ~ 5> fc z- 2fe , Q{z) ~Y.1kZ- 2k -\ 

k=0 fe=0 
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with some (explicit) coefficients Pk, Qk- In particular, p\ — 1 , qi — —1/8. Then 
(2tt) 2 J 2 (z) ~ — ((P + \Q)S"-*M + (P - iQ)e-^-^Y 

~ — (-i(P + iQ) 2 e 2i * + i(P - iQ) 2 e- 2iz + 2P 2 + 2Q 2 ) . 

We set 

(fO.lf) A = AttID^ 1 {P 2 (\D\) + Q 2 (\D\)) , F ± = T 27ri\ D \- 1 (P(\D\) ± iQ(\D\)) 2 e ±2i ^ . 

This completes the proof. □ 

We will now connect this to Theorem |2.f| Let p a — (0,0), go = (2,0), v — (0,7r), wo — (0, 7r). Then 
v a e S(p ). Choose £ = (1,0), conormal to the conjugate locus E(g ) = {\x — q \ = 2} at p ; and choose 
i]o = (1, 0), conormal to the conjugate locus S(p ) = {\x — Po\ — 2} at q . The directions of £o, Vo reflects 
the choice of the orientation we made earlier. We refer to Figure [5] 




Figure 5. 



If we localize R near v — vq, then the pseudo-differential part of R*xR is (1/2)A), see (5.10). Therefore, 



in the notation of Theorem 2.1 



.4 



1 



'-Ac, F = F + 



The canonical relation of F + maps (po,£,o) m t° (lOiVo)) see Figure [5j while that of P_ maps (po,—£,o) into 
(qo,~Vo)- This is consistent with the results in Theorem 2.1 where the Lagrangian has two disconnected 
components located near (po, go, ±£o, T^o)- 

To analyze the operator ( |9.9| , note first that A x = A 2 = A /2. Let us first analyze this operator applied 
to distributions with wave front set near (go, rjo) but not near (go, — f]o)- Then F reduces to F + only, and we 
have, modulo vf -00 , 



AJ 1 F*A^ 1 F = ^A- 2 F^ 



Id, 



see (10.11). The analysis near (go, — tjq) is similar. Therefore, we have a stronger version of Theorem 9.2 
this case: singularities can cancel to any order. 
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Theorem 10.2. Let fx be any distribution with WF(/x) supported in a small conic neighborhood of some 
(xq, £ ) G T*R 2 \0. Then there exists a distribution /2 with WF(/2) supported in a small conic neighborhood 
of (xq ± 2£°/|£°|, that is an image of WF(/i) under the map T±, so that R(fi + f%) G C°° for all unit 
circles in a neighborhood of the unit circle C(xq ± £ ). 

In other words, for a fixed circle Cq of radius 1, there is a rich set of distributions /, with any order of 
singularity at AT* Co, so that those singularities are invisible by X localized near Co, i.e., Xf G C°°. Explicit 
examples can be co nstru cted by choosing /^(s) = S(x — qo), then F fi near po is just given by the Schwartz 



kernel of R*R, see (10.6). To obtain f±, we apply 2Aq 1 to the result. 

We would like to emphasize on the fact that the theorem provides an example of cancellation of singularities 
for the localized transform only. As we will see below, Rf € C°° (globally) for f E £' implies / G C°°. 
On the other hand, without the compact support assumption, on can construct singular distributions in the 
kernel of R, using the Fourier transform. 

10.1.4. The wave front set of a distribution in Keri?. Now, if Rf = or more generally, if Rf € C°°, one 
easily gets that 

(10.12) V/ G Keri?, WF(/) is invariant under the action of the group {J 7 ™, m £ Z}. 

Then, if / is compactly supported (or more generally, smooth outside some compact set), we get that WF(/) 
must be empty, i.e., / G C°°(R 2 ). In other words, even though recovery of WF(/) is impossible by knowing 
Xf locally, as we saw above; the condition Xf G C°° globally, together with the compact support assumption 
yielded a global recovery of singularities. Here an important role is played by the fact that X is translation 
invariant, and in particular, our assumptions are valid for any (po-,&o) G TSTi, 2 that cannot be guaranteed 
in the general case. Also, the dynamics is not time reversible; therefore for each (xq,^°) G T*M\ there 
are two different curves through xq in our family. The latter is true for general magnetic systems with a 
non-zero magnetic field, see [3]. 



Remark 10.1. One can see that R is invertible on L 2 (M) by using Fourier transform, see (10.2 1. The 
formal inverse is 1/J (|£|), and conjugating a compactly supported \ with the Fourier transform, one gets a 
convolution in the £ variable that will smoothen out the zeros of Jo(|£|), thus producing a Fourier multiplier 
with asymptotic ~ ICI 1 ^ 2 - 111 L P (R 2 ) with p > 4 however it is not invertible, and elements of the kernel 
include functions with Fourier transforms supported on the circles Jo(ICI) — 0> see a l so [T^l ITj. 

Finally, we remark that in this case, one can study R directly, instead of R*R = i? 2 , with the same 
methods. Our goal however is to connect the analysis of this transform with our general results. 

10.2. The X-ray transform on the sphere. Consider the geodesic ray transform on the sphere S n . The 
conjugate points are not of fold type, instead they are of blow-down type. Let J be the antipodal map. 



Without going into details, we will just mention that then (2.3) still holds with 



CN = ID]- 1 - ID]- 1 J, 

with some constant C, where the canonical relation of F is the graph of the antipodal map, lifted to T* S 2 . 
Then CiV|D| = Id — J. The canonical graph is an involution, however (its square is identity), so arguments 
similar to that in the previous example do not apply. That means that singularities may cancel. In fact, it 
is known that R has an infinitely dimensional kernel — all odd functions with respect to J. This is a case 
opposite to the one above. 

In this case S consists of all antipodal pairs (x,y), and has dimension 2 (and codimension 2), unlike the 
case above (dimension 3 and codimension 1). On the other hand, A/"*S still has the same dimension (that is 
2n=4, and this is always the case as long as S is smooth submanifold) . One can see that the Lagrangian in 
this case is still Af*Y,. 
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10.3. Magnetic geodesies in R 3 . Consider the magnetic geodesic system in the Euclidean space R 3 with 
a constant magnetic potential (0, 0, a), a > 0. The geodesic equation is then given by 

(10.13) 7 = 7 X (0,0, a), 

where x denotes the vector product in R 3 . The r.h.s. above is the Lorentz force that is always normal to 
the trajectory and in particular does not affect the speed. We restrict the trajectories on the energy level 1 
that is preserved under the flow. Then we get 

7 1 = a^ 2 , j 2 — — aj 1 , j 3 = 0. 

The magnetic geodesies are then given by 

7(i) = 7(0) + (-(sm(at + 9) - sin6>), -(- cos(at + 6) + cos0), tz) , 

where (r, 9, z) are the cylindrical coordinates of 7(0). The unit speed requirement means that 

r 2 + z 2 = 1. 

The geodesies are then spirals; when z = then they reduce to closed circles, and when r = they are 
vertical lines. 

The parameterization by cylindrical coordinates is singular when r = 0. Away from that we can use 9, z 
to parametrize unit speeds. Then in exp p (y), we use the coordinates (t, 9, z) to parametrize v, i.e., 



tU/l - z 2 (cos6»,sin6>),. 



At t — we may have additional singularity but this is irrelevant for our analysis since we know that the 
exponential map has an injective differential near v = 0. An easy computation yields that the conjugate 
locus is given by the condition at — 7r, i.e., 

S p (v) = {v; \v\ = ir/a} , 

and this is true for any p € R 3 . This is a sphere in TR™. For S p we then get 

(10.14) 7(7r/a) = p + a~ 1 (-2r sin6, 2r cos6», ttz) 

with p = 7(0). This shows that is an ellipsoid 

' -2 1 1 /„ „ \2 , 1 i„ „ \2 -2 



s = i(P,Q)\ j (9l -Pi) +|fe-P2) +^(93-^3) =a 
Then 

(10.15) AA*E= |(p,g^,77); ( M )eS; g = c(p x - q u p 2 - q 2 , ^( P3 - q 3 j) , V = ~t, c e r| . 

Therefore, given p 7 we can immediately get q as a smooth function of (p, £), and we can obtain w so 
that exp p (w) = q by (10.14), where the l.h.s. is g. Therefore, (p, £) 1— > u is a smooth map, and therefore 
(P>£) ^ ( ( ?7 ? 7) i s a smooth map, too. The later also directly follows from (10.15), since 77 = — £. 
We therefore get that F is an FIO of order —3/2 with a canonical relation 

(10.16) (p,0->(9»0, 

where g can be determined as described above. A geometric description of q is the following: q is one of the 
two points on the ellipsoid S, where the normal is given by £. The choice of one out of the two points is 
determined by the ch oice of the initial velocity v near which we localize; changing v to —v would alter that 
choice. Since (10.16) is a diffeomorphism, F is of canonical graph type, and therefore maps H s to fp+ 3 / 2 . 
In contrast, A\^, are elliptic of order — 1, thus they dominate over F. By Corollary 9.1 X can be inverted 
microlocally in the setup described in Section [2j 



30 



P. STEFANOV AND G. UHLMANN 



10.4. Fold caustics on product manifolds. Let (M,g) = (M',g') x [M" ,g") be a product of two Rie- 
mannian manifolds. The geodesies on M then have the form 

Consequently, 

expp(w) = (expp,(u'),exp£„ (?;")). 

Assume that in (M',g'), v is conjugate at po of fold type, and assume that v' ' is not conjugate at p ' in 
(M",g"). Then 

dexp p (w) = diag(dexp p ,(t/),dexp p „(i/')). 
The kernel of dexp p (w) then consists of N p (v) = N p ,(v') x 0. Next, S(p) = S(p') x T p »M", and £(p) = 
E'(p') x M" . Then N p (v ) is transversal to S(p) at v = Vq, therefore (v',v") is a fold conjugate vector 
for v' £ S'(p) close to vo and for any v" . Then the left projection 7Tl of the Lagrangian J\f*Jl consists of 
(p,£) with S 7r L (E') and £" = 0. Thus the rank drops at least by n" = dim(M"). We get the same 

conclusion for 7Tr(.A/"*I!). Therefore, A/"*S is not a canonical graph in this case. 

Let n' = dim(M') = 2. Then the canonical relation in (M',g') is a canonical graph, and we get that 
7Tl have rank 2n' + n" = 4 + n" instead of the maximal possible 2n = 4 + 2n"; i.e., the loss is exactly 

n"'. 

Assume now that nf = 2, n" = 1, and the metric in M is given by 

2 

5^ ^(^a^dx'Ma;' 3 + (dx 3 ) 2 . 

Q,/3=l 

Assume also that in M', we have a fold conjugate vector vq — (0,1) at x 1 — x 2 — 0. Then all possible 
conormals to the conjugate loci at (0, 0) corresponding to small perturbations of v will lie in the plane 



v 3 = 0. This is an example where Corollary 9.2 can be applied. We can recover singularities of the kind 
£ = (^1,^2,^3) at po = (0,0,0) with £ 3 ^ and (£i,£a) m a conic neighborhood of (1,0). The ones with 
£3 = are the problematic ones. 
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